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Abstract 

The dynamics of a molecular junction consisting of a PTCDA molecule between the tip of a 
scanning tunneling microsope and a Ag(lll) surface have been investigated experimentally and 
theoretically. Repeated switching of a PTCDA molecule between two conductance states is studied 
by low-temperature scanning tunneling microscopy for the first time, and is found to be dependent 
on the tip-substrate distance and the applied bias. Using a minimal model Hamiltonian approach 
combined with density-functional calculations, the switching is shown to be related to the scattering 
of electrons tunneling through the junction, which progressively excite the relevant chemical bond. 
Depending on the direction in which the molecule switches, different molecular orbitals are shown 
to dominate the transport and thus the vibrational heating process. This in turn can dramatically 
affect the switching rate, leading to non-monotonic behavior with respect to bias under certain 
conditions. In this work, rather than simply assuming a constant density of states as in previous 
works, it was modeled by Lorentzians. This allows for the successful description of this non- 
monotonic behavior of the switching rate, thus demonstrating the importance of modeling the 
density of states realistically. 
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I. INTRODUCTION 



The scanning tunneling microscope (STM) is a valuable and versatile tool for the study 
and manipulation of nanoscale structures.- 1 ^ In scanning mode, it can be used to image 
surfaces with atomic resolution, and to probe the electronic density of states at a range of 
energy values. Alternatively, it can be brought into contact with surface features to form 
junctions and measure transport properties.- - - Nanostructures and devices can be manip- 
ulated and fabricated using an STM, with the possibility to pick up and deposit atoms 
and molecules using the tip.-^~— An important aspect related to the tip-molecule inter- 
action is the telegraph noise observed in the conductance in certain circumstances, which 
originates from the repeated switching of single atoms or functional groups between differ- 
ent stable configurations.-^— Several physical mechanisms have been proposed to explain 
this phenomenon: thermal activation, vibrational heating (for intermediate biases yi2r21 and 
transition through an electronic excited state with no conformational bi-stability (for high 
biases).— If the masses involved are not too large (i.e. for a single atom), quantum tunneling 
is also possible.— 

In this work, we present a systematic study of this switching behavior in the specific system 
of perylene-3,4,9,10-tetracarboxylic-dianhydride (PTCDA, inset in Fig. QJa)) on Ag(lll), 
using both experimental and theoretical methods. With results from density- functional 
(DFT) calculations and by extending a microscopic model developed in Ref.Qto describe 
the coupling of an adsorbate energy level to the adsorbate vibrational excitations, a good 
agreement with the experimentally measured switching rates can be achieved. 

PTCDA deposited on Ag(lll) forms a highly ordered metal-organic interface, the elec- 
tronic and geometric structure of which has been well-characterized using a variety of both 
experimental and theoretical techniques.— ~— The PTCDA molecules form long-range or- 
dered commensurate monolayers on the Ag(lll) substrate with two flat-lying chemisorbed 



molecules per unit cell in a herringbone arrangement (see Ref. l25l ). The chemisorption re- 
sults in the former lowest unoccupied molecular orbital (LUMO) of the isolated molecule 
being shifted below the Fermi level of the silver surface, so that there is charge transfer from 
the substrate to the molecule thus producing a net negative charge on the molecule.— 

In previous experiments we have found that it is possible to form a chemical bond between 
the carboxylic oxygen atoms and the STM tip, if the latter is approached towards the 
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FIG. 1. (color online) (a) Schematic picture illustrating the up- and down-states and the switching 
between the two for the tip-PTCDA-Ag(lll) junction (red-oxygen, yellow-carbon, light-blue- 
hydrogen, grey-silver). The atomic coordinates are taken from the DFT calculations described in 
Ref. [g. The inset shows the structure of the gas phase molecule, (b) Measured current at 125mV 
during approach of the STM tip above the carboxylic oxygen of the PTCDA molecule in (c). The 
tip was moved by 0.6A at a rate of lA per 23 min. (c) STM image of the edge of a monolayer 
of PTCDA. The white arrow indicates the PTCDA molecule which was used for the switching 
measurements, and points to the oxygen atom which interacts with the tip. 



molecule above one of the carboxylic oxygen atoms:— >^ the oxygen atom, followed by part 
of the carbon skeleton of the PTCDA molecule, jumps into contact with the tip. The most 
likely distance for this single switch to happen (without applying a bias voltage) is 6.65A.- 
In a theoretical analysis, carried out by calculating potential profiles of relaxed PTCDA 
molecules between tip and surface as a function of oxygen-surface separation for a range tip 
sample separations, we found the spontaneous jump into contact at 6.2A,- in good agreement 
with experiment. 

Once the molecular junction with the tip has been formed, there are two possible ways for 
the molecule to behave when the tip is retracted: either the molecule is peeled off from the 
surface completely or it falls back to the surface.- We have further observed that, under 
certain conditions (see below), the current fluctuates in time between a high- and a low- 
conductance state, see e.g. Fig. QJb) in which the telegraph noise in the current is evident. 
These two-state fluctuations can be explained by the switching of the molecule in and out 
of contact with the tip (see Fig. Ufa)). In the high- conductance state, one of the carboxylic 
oxygen atoms of the molecule forms a chemical bond with the tip ("up-state"), establishing a 
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two-terminal molecular junction, while in the low-conductance state the molecule is bonded 
exclusively to the surface ("down-state") so that a tunnel barrier is now present between 
the tip and the molecule. These switching processes of the molecule can also be seen in the 
topographic images taken with the tip very close to the surface.- 



II. EXPERIMENTAL METHODS 

Our experiments have been performed with a CREATEC low temperature scanning tunnel- 
ing microscope (5-6 K) in ultrahigh vacuum with a base pressure below 10~ 10 mbar. The 
Ag(lll) surface has been prepared by repeated sputtering/annealing cycles (Ar + ion energy 
0.8 keV, annealing at approximately 850 K). Surface quality has been controlled in situ with 
low-energy electron diffraction (LEED). The PTCDA molecules have been evaporated from 
a Knudsen cell at 580 K onto the surface at room temperature. An electrochemically etched 
tungsten wire has been used as the STM tip, which has been cleaned in situ by annealing. 
The final atomic sharpening has been done by the indentation of the tip into the clean 
metal substrate and/or by the application of voltage pulses. Tip quality has been checked 
by measuring the surface state of Ag(lll). The PTCDA material (commercial purity 99%) 
has been purified by resublimation and outgassing in ultra high vacuum. 
Prior to the measurement of the switching process, the STM tip was stabilized at Vbi as =- 
340mV and 1=0. InA, corresponding to a tip-surface separation of 10. 6A (Ref. which 
is outside the regime in which repeated switching is observed. Absolute calibration of the 
tip-surface separation was done as described in Ref. |8| (error of ±0.5A for the absolute 
height). Time spectra of the current were recorded for different bias voltages and tip-surface 
separations with the feed-back loop switched off. The time dependent current I(t) is shown 
in Fig. [2](a) for the applied bias voltage of 95mV and with the tip positioned at 7.1 A above 
the substrate. 

The quantitative analysis of the switching process, which is the primary objective of this 
paper, has been carried out for molecules located at the edge of a monolayer island of 
PTCDA/Ag(lll) (as indicated with the white arrow in in Fig.[T^b)). The reason for choosing 
these molecules is that the PTCDA molecules in the midst of a compact layer are more 
difficult to pick up due to strong intermolecular interactions with neighboring molecules via 
hydrogen bonds, 2 - while isolated molecules do not always fall back to the same position on 
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FIG. 2. (color online) Switching of a PTCDA molecule between up- and down- state, (a) Current 
vs. time trace measured at 95mV with a tip-surface separation of 7.1 A. (b) Map of the average 
switching frequency as function of bias voltage and tip-surface separation. The corresponding 
spectra were measured at constant bias during tip approach. The bias range from -120mV to 
120mV was covered with a step of 5mV. (c), (d) Residence time histograms for the up-state and 
the down-state, extracted from the time trace in panel (a). The red solid lines show the exponential 
fit used to extract the transfer rate R. 

the surface when they switch from the up- to the down-state, thereby leaving the junction 
and precluding the continued measurement of the switching time trace. 
A color-coded map of the frequency of switching events as a function of bias voltage and tip- 
surface separation is displayed in Fig.[2](b). We observe the following: (1) Repeated switching 
occurs for both bias voltage polarities above a threshold of approximately |100|meV. In con- 
trast, for Ubias < |100|meV a single jump into contact occurs^ (not indicated in Fig. 12(b)); 
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for these bias voltages, the junction may only (but does not necessarily) switch back from 
the up- to the down-state if the tip is retracted again beyond the tip-surface separation at 
which the jump into contact has originally occured (hysteresis). (2) Repeated switching 
occurs in a narrow bracket of tip-surface separations in the range from 7.34A to 7.14A. 
(3) The range in which repeated switching is observed appears at slightly larger tip-surface 
separations for negative bias than for positive bias. This latter fact may be related to the 
negative polarization of the carboxylic oxgen atoms in Ag(lll)-adsorbed PTCDA. 
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FIG. 3. (color online) Double logarithmic plot of the transfer rate for PTCDA switching between 
the STM tip and a Ag(lll) surface for different tip-surface separations. Measured transfer rates 

symbols. Solid 



22. The inset in 



for switches from (a) surface to tip and (b) tip to surface are indicated by smal 
lines display the theoretical transfer rate, (a) dashed orange lines the model of Ref. 
(b) shows a possible fit with the model of Ref. 22. However, the parameters thus obtained disagree 
with both experiment and DFT calculations. 



iFrom the I(t) curve in Fig. |2ja) one can see that for the chosen bias and tip-surface 
separation, the up-state is preferred: the statistical residence time analysis (Figs. [2^c) and 
(d)) reveals a difference of more than one order of magnitude in the residence time values 
for the high- and low- conductance states. The single exponential behavior of the curves 
indicates a two state Markovian switching process where the residence time probability 
density P is given by the expression P (t) = R exp (— Rt) . Here, R is the transfer rate 
between the two conductance states. It is obtained by fitting the equation for P to the 
corresponding residence time histogram. By performing such a transfer rate analysis for 

7 



different bias voltages one can determine the transfer rates as a function of bias for every 
measured tip-surface separation. The rates for three typical tip-surface separations are 
displayed in Fig. [3J The tip— ^-surface transfer rate increases monotonically with applied bias 
in the given voltage range, but the surface— >tip transfer rate appears to have a maximum 
around 180mV. Finally, both rates are dependent on the tip-surface separation, as can also 
be seen in the experimental data in Fig. |2^b). 

III. THEORETICAL MODEL 

To gain insight into the observed current switching, we first focus on the nature of the 
coupling between PTCDA and the surface, and then provide a link to the experimental data 
by applying a model calculation. The mechanisms of the chemical bonding of PTCDA to 
Ag(lll) includes hybridization of the molecular orbitals with the substrate states, charge 
transfer between the substrate and the molecule, local bonds of the carboxylic oxygens to 
silver atoms below and an extended bond of the molecular 7r-system to the surface.- 1 ^— 
Assuming that the two meta-stable positions can be well-represented by a (not necessarily 
symmetric) double-well potential, the transfer of an adsorbate between the two minima 
may involve a variety of physical processes, such as (i) thermal activation, (ii) quantum 
tunneling, (iii) a transition through an electronic excited state with no conformational bi- 
stability, or (iv) vibrational heating. Process (i) is of minor interest in this work, since the 
experiments are performed at very low temperatures (5-6 K) and the barrier height is larger 
than lOOmeV, which excludes the thermal activation. Due to the relatively large mass of the 
part of the molecule involved in the switching process, process (ii) is also very improbable. 
Assuming a tunneling barrier of lOOmeV height (measured from the vibrational ground 
state) and 1A width (cf. Fig. H^a)) the corresponding tunneling rate for the carboxylic 
oxygen atom was estimated to be of the order of 10 _8 Hz. For process (iii), which involves 
an excited state of the molecule, the residence time of the tunneling electrons has to be 
sufficiently large to induce this excitation. However, since the molecule is chemisorbed on 
the Ag(lll) surface, this residence time is expected to be quite small, so that process (iii) 
also seems unlikely in this case. Thus, we suggest that the microscopic mechanism leading 
to switching is related to vibrational heating, where the transition is induced by progressive 
vibrational excitation of the relevant chemical bond (i.e., either the oxygen-surface bond 
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FIG. 4. (color online) (a) Double well potentials of relaxed PTCDA molecules between tip and 
surface as a function of oxygen-surface separation for a range tip-sample separations calculated with 
DFT for a single PTCDA molecule.— (b) Schematic double well potential used to describe up- and 
down-states of PTCDA in the STM junction, including vibrational levels and model parameters 
appearing in Eq. ([6]). IY^- denote the relaxation and excitation rate of a molecular vibration due 
to scattering of tunneling electrons, and n\ and n-i denote the critical number of vibrations which 
have to be excited to induce the switching. 



for the surface— >-tip process or the oxygen-tip bond for the reverse process) by the inelastic 
scattering of tunneling electrons, eventually leading to bond breaking. The transition rate is 
then mainly determined by the competition between energy gain from the tunneling charges 
and energy losses due to electron- hole pair generation and/or coupling to the substrate 
phonon continuum. 

In Ref. |22| Gao et al. developed a theoretical model to describe atomic switching by 
vibrational heating. They concluded that the switching rate should exhibit a power law 
dependence R oc Vuas 1 on the bias voltage Vuas where n is the number of vibrational levels 
that have to be climbed before the switch can occur. In our experiments we observe a 
striking difference between the tip— ^-surface and surface— )-tip switching processes as far as 
the bias dependence of the switching rate is concerned (cf. Fig. |3]). While the tip— ^surface 
process shows an almost linear RjVuas) behavior in the double logarithmic plot, in essential 



agreement with the prediction of Ref. |22|, a reasonable description within the model of Gao 
et al— for the surface— >tip process is very unlikely, because RiVuas) deviates from a simple 



power law, showing a saturation of the transfer rate at approximately 120mV, with even 
a possible maximum around 180mV. Moreover, a (rather poor) fit of the data in Fig. [3(a) 
with R oc Vbi as n would only be possible by assuming n = 1 (cf. dashed line in Fig. [3(a)), 
which is in contrast to the calculated potential energy surfaces, see Fig. 11(a). 
Below we show that if the energy dependence of the density of states around the Fermi 
level is taken into account explicitly, and if in particular different transport orbitals for the 
two configurations (i.e. up- and down-states) are used, the evident differences between the 
surface— >tip and tip— ^-surface processes can be rationalized and both the data in Figs. [3(a) 
and (b) can be fitted with parameters which are in qualitative agreement with DFT results 
for the electronic structure of the molecular junction. In our model, we will neglect the 
coupling to the substrate phonon continuum, since the anharmonic coupling is, in general, 
very small at low temperatures. 

Our model is a minimal approach based on that used in Ref. 
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to describe the vibrational 



heating. The Hamiltonian describing the tunneling of electrons between the STM tip and 
the surface via an adsorbate level (in this case, the adsorbate being the PTCDA molecule) 
has the following form: 

H = } j £ s c\c s + E t c\c t + EmC^Cm + hu b ] b (1) 

s i 

+ 

s i 

Here s, t and m label one-electron states \s), \t) and \m) of the surface, the tip and the 
molecule, respectively, with the corresponding energies e s , e t and e m . The hopping between 
the surface and the tip via the molecular level is described by the two terms including T sm 
and T tm . The coupling between the vibrational motion of the molecule and the electron 
propagating through it can be modeled by: 

H e _ v =\ (&t + b) (4c m ) , (2) 

where, 



A ° = V2iL^ ■ (3) 

The coupling is modeled by assuming that e m is a linear function of the vibrational coordinate 
q, £ m {q)'-, u) is the frequency of the molecular vibration with the normal coordinate q = 
\J (2Mw) + ^) an d mass M, and e' m = de m /dq at q = 0. 
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Since the effect of the electron- vibration interaction on the adsorbate electronic states is in 
general weak, it can be treated by first-order perturbation theory. The assumed linearity in 
the charge-vibron coupling simplifies the problem since only the excitation and relaxation 
rates, and T^, between the vibrational ground state and the first excited state are required 
(cf. Fig. IHb)). In first-order perturbation theory these transition rates are given by Fermi's 
Golden Rule: 

r t =2 Y ElO'' 1 !^-^' )! 2 fi^-fj) Siej-et + tiu), (4) 
T ± =2 T Z}lfr°l#e-«|U)| 2 Ml-/;) 5{e j -e l -f^) (5) 

3,1 

where and 1 are the vibrational ground state and the first excited state respectively, while 
j and I denote any of the stationary one-electron states of the tip or the substrate with corre- 
sponding Fermi-Dirac distributions f^i = 1/ {1 + exp [(e — £jj)/(/cbT)]}, and H e - V denotes 
the electron- vibration interaction (Eq. (j2J)). 

These rates describe the vibrational excitation and relaxation induced by the tunneling 
electrons. Since the initial and final states of a tunneling electron can be located either in 
the tip or the substrate, these rates can be decomposed into four different terms: Fi s , , rl* , , 
r|*|_ and T^, which sum up to give T^. Here, the first (second) superscript denotes whether 
the final (initial) state belongs to the surface or the tip. In contrast to Ref. [22j, we will not 
assume that the adsorbate local DOS is constant over the relevant energy range, but rather 
we model it by a Lorentzian shape, p s ^(E) = A Sit /((E — e m ) 2 + A 2 ), where A = A s + A t , 
with A s and At describing the coupling between the molecular level and the substrate and 
tip electronic states, respectively. Using this function, the excitation and relaxation rates 
can be calculated analytically in the low temperature limit. We refer the interested reader 
to appendix A for further details and a comprehensive description of the calculation. 
To describe the transfer between the two possible meta-stable states a truncated harmonic 



oscillator model, as described in Ref. 



221 . is adopted. The transfer rate R can be expressed as 



a product of the transition into level n (see Fig.|U(b)) and an effective Boltzmann factor (with 
characteristic temperature T v = fiw/ (kB^lPi/T-f])) describing the probability to arrive at 
the sub-critical level n — 1 where the transition takes place:— 



R ~ n r+ exp 



(n — 1) hid 
k B T v 



nr t =L . (6) 
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Since the adsorbate local DOS is not assumed to be constant over the relevant energy range, 
the above expression in general does not yield a simple power law dependence on the applied 
(R oc Vbi as n ). This simple scaling law can only be recovered, if the molecular 



bias as in Ref. 



22 



level is situated far from the Fermi energy (so that the DOS at Ep is almost constant) 



IV. RESULTS AND DISCUSSION 

Using Eqs. (A1-A4) in the appendix, we are now able to fit the transfer rate in Eq. ([6]) to the 
experimental results. Figure |3] shows the fitted transfer rates as a function of bias voltage, 
together with the experimental data. The corresponding fitting parameters are listed in 
Tab. [J and will now be discussed in detail. 

The vibrational energies Hou (i.e. the size of the steps on the "vibrational ladder") were 
determined from the curvature of the calculated potential energy surfaces, shown in Fig.0Ja). 
They lie around 19meV for the shallower well of the down-state, and around 40meV for the 
deeper well of the up-state. The dependence of these vibrational frequencies on the tip- 
surface distance is negligible (cf. Tab.|T]). 

The n are an output of the fitting of the transfer rates. Multiplied with tkj, they yield 
the barrier heights for the switching process. The products riihui and ^2^2 in Tab. U 
are consistent with the potential energy surfaces obtained from DFT calculations shown 
in Fig. HJa), which exhibit a highly asymmetric double well, with a shallow well for the 
down-state and a deep one for the up-state. The asymmetry increases as the tip-surface 
separation is reduced. In particular, the depth of the potential well of the up-state (nifuvi), 
which according to Tab. [J amounts to 0.53eV at 7.3A, agrees quite well with that calculated 
within DFT, whereas the model predicts a down-state well of 0.17eV at 7.3A that is slightly 
deeper than that derived from the ab-initio calculations (cf. Fig. H(a)). This may be due 
to the fact that the potentials in Fig. 0Ja) were calculated for a single PTCDA molecule, 
whereas in the switching experiments edge molecules were used; their hydrogen bonds to 
neighboring molecules will lead to a significant increase of the barrier height. Note, however, 
that the model does correctly predict the decrease in depth of the down-state well as the 
tip-surface separation is decreased; this tendency is due to the reduction of the potential 
minimum to a saddle point for tip-surface separations of less than about 6.2A (cf. Fig.|4](a)). 
A further important parameter in our model for the transfer rate is the position of energy 
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level e m through which the electron current that causes the vibrational heating passes (i.e. the 
transport level), because this influences the energy dependent density of states that enters 
the rate via Eqs. (jl]) and (jSJ). It is clear that levels on either side and closest to the Fermi 
energy Ep are the most important channels for the electron current. Our DFT calculations 8 
show that mainly states both above and below the Fermi level could contribute, see Fig. 
The level below Bp is the former LUMO that gets filled on adsorption and that is clearly 
observed in scanning tunneling spectroscopy— ~— The sharp level above Ef that is found in 
DFT appears in experiments as a broader feature in the gap between the former LUMO and 
the LUMO+1, especially for molecules at the edges of monolayer islands. In our minimal 
model Hamiltonian we can only take one transport level into account. It turns out that the 
qualitatively different behavior of the two processes surface— )-tip and tip— ^surface requires 
the use of two different transport levels, depending on the switching direction. This is 
reflected in Tab.Uby negative values £ m> i for the tip— ^-surface process, while the surface— ?-tip 
process has positive e mj2 values (the spectral density of the levels e m ^ and £ m>2 are shown 
in Fig. [6]). In other words, we have to assume that in the up-state the switching current 
passes mainly through occupied DOS of the junction, whereas in the down-state it passes 
predominantly through the empty DOS of the adsorbed molecule. Note that due to the way 
in which the bias voltage drops between tip and substrate, both molecular levels e m) i and 
£m,2 are within the bias window and may in principle contribute to the transport, but in our 
minimal model we can - as mentioned above - only take one into account at a time. 
The fitted values £ TO) i show a clear tendency to move up towards the Fermi level as the tip- 
surface separation is increased. This tendency is known very well both from experiment^ and 
DFT calculations,-^ although the precise level positions in experiment and ab initio theory 
differ from those in Tab. [H This is not too surprising since our minimal model only allows 
for a single Lorentzian level whereas the actual density of states is much more complicated. 
The fitted values e m ^ range between 0.24eV and 0.30eV, whereas the DFT calculation has 
this level fixed at 0.2eV. 

The small decrease in the transfer rate of the surface— s-tip process at about 220meV 
(Fig.[3]^a)) is due to the molecular level e m ^ entering into resonance with the Fermi energy of 
the STM tip, which leads to a reduction of the vibrational lifetimes of the PTCDA molecule 
in the junction {i.e. the rate Y^ at which the molecular vibrational energy dissipates into 
the electrodes is increased, cf. Eqs. (Al, A4)). This in turn reduces the transfer rate of the 

13 




for the PTCDA molecule on the Ag(lll) surface (down-state-dashed red) and attached to the tip 
(up-state-black) for a tip-surface separation of 7A. The level just above the Fermi energy is at the 



same position as in the simple model described here. The level below, however, is lower in energy 
compared to the model but also compared to the experiments. 



T 1 1 1 r 




(E-Ep) [eV] 

FIG. 6. (color online) Density of states obtained from the fitting procedure for different tip-surface 
separations. If the molecule is attached to the tip the level below the Fermi energy moves up with 
increasing distances while the level above only shifts slightly. 



14 



molecule between the surface and the tip. Further raising the bias voltage beyond this point 
results in the transfer rate increasing once again due to non-resonant tunneling. We stress 
that this behavior can only be obtained if an energy-dependent DOS is used; a constant 
DOS could not yield such a behavior. Unfortunately, the increase above 240meV cannot be 
observed in the experiments since the molecule normally disintegrates at lower biases than 
this because of the high current density. 

V. CONCLUSION 

In summary, switching between low- and high conductance states has been observed in 
a single molecule junction consisting of a PTCDA molecule on a Ag(lll) substrate and 
contacted by an STM tip. The rates for the transition between these two states can be 
sensitively tuned by varying the applied bias as well as the tip-surface separation. A vi- 
brational heating mechanism where molecular bonds are excited by tunneling charges has 
been proposed to interpret the experimental results. Switching rates were calculated within 
a minimal model Hamiltonian approach describing the interaction between tunneling elec- 
trons and local molecular vibrations. The experimental results could be fitted over a broad 
voltage range for the cases where the PTCDA molecule switches both from the surface to 
the tip and from the tip to the surface. In particular, the non-monotonic behavior of the 
surface to tip switching rate could only be described by modeling the DOS by Lorentzian 
functions instead of assuming it to be energy independent, as has been the usual practice 
in the literature until now. This demonstrates that it is crucial to take the non-constant 
behavior of the molecular DOS into account. 
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TABLE I. Model-parameters for the switching of the PTCDA between the tip to the surface 
obtained by fitting R to the experiments. The subscripts "1" and "2" indicate switching from 
STM tip to the surface or the reverse process, respectively. Energies e m , Hoj, nhoj, A s and A t 
are given in meV. The dimensionless parameters Ai,2 = Ao/^1,2 describe the electron- vibration 
interaction. 
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0.025 


-136 10 40.94 


409 


148 


33.59 


7.44A 


0.024 


-124 


9 


40.98 


369 


148 


26.30 


7.47A 


0.010 


-110 


8 


41.00 
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148 


26.87 


Tip-surface distance 


A 2 


£m,2 


n 2 


Huj 2 




A s , 2 


A t , 2 
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260 


8 


18.95 


152 


24 


11.9 


7.20A 


0.009 


257 


8 


19.00 


152 


25 


12.5 


7.24A 


0.007 


253 


9 


19.06 


172 


30 


14.9 


7.27A 


0.006 
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9 


19.12 


172 


28 


13.9 


7.30A 


0.005 


260 


9 


19.15 


172 


24 


11.8 


7.34A 


0.004 


259 


9 


19.22 


173 


22 


10.8 


7.37A 


0.003 


269 


10 19.26 


193 


23 


11.4 


7.41A 


0.002 


285 


11 


19.32 


213 


23 


10.1 


7.44A 


0.002 


299 


10 19.37 


194 


10 


4.9 


7.47A 


0.0004 


258 


10 19.42 


194 


28 


8.9 
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Appendix A: Calculation of the transition rates 



In the following we want to sketch the derivation of the transition and relaxation rates. The 
terms and , are all similar, and it is sufficient to calculate explicitly only the term 
TJ S . Inserting the electron- vibration interaction (Eq. ([2])) into Eq. ([5]) together with the 
expression for the molecular DOS gives 

-/ ^2 



7T £ 
-pss 2 v " "' 



Mu 

Mu 
A 2 (e' ) 2 

o s \ mJ 



^ \{a'\m)(m\a)\ 2 [l - f s (e a >)] f s (e a ) 5{e a > - e a - tuu) 



Mu 



7T 



A 2 ie' ) 2 

M U 7T 



dB^(e)AC(e + /iw) [l-/.(e + fiw)l /.(e) 



[e-e m ] 2 + A 2 [e + ^-e m ] 2 + A 2 



de 



[l-f s (e + hu)} f s (e) 



e - E, r 



A 2 [e + fru- e m Y + A 2 



x [1 — Q(e Fs — Hu — s)] Q(e Fs — e) 



A 2 is' ) 2 
2 8 K m) I de 

Muix 

Epg—huj 



l 



e - £mf + A 2 [e + fku- e m } 2 



A 2 



4 A 2 A 2 



7rAfr 2 w (4 A 2 + ^ 2 o; 2 



hu ( tan 1 



e Fs + hu 



A 



tan 



£m — £Fs — hu 

A 



+ A (log [A 2 + (e m - £ F s + ^) 2 ] 



+ log [A 2 + (e m - e Fs - hu) 2 } - 2 log [A 2 + {e m - e Fs ) 2 } ) 



(Al) 



In the first step the sum over states has been replaced with an integral over e by introducing 
p s m {e). In the second step the expression for the molecular DOS was used to rewrite the local 
density of states. Since the STM experiments are carried out at 5-6 K one can approximate 
the Fermi function with the Heaviside step function in the next step. Thus, the limits 
of the integral can be changed from +oo and — oo to e Fs or e Fs — hu respectively. The 
influence of an applied bias can be easily introduced by shifting the Fermi level of the 
surface e Fs = e F o s + eV, where e F Q S is the Fermi level at V = of the surface. Since we used 
the low temperature approximation in step 3 of Eq. (Al) the excitation rates Ti s ' tt become 
zero, because of the Pauli exclusion principle. The Pauli exclusion principle also simplifies 
the calculation of the remaining terms and ri*, , which describe the transition rates 
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due to the inelastic scattering of tunneling electrons between surface and tip. Assuming 
£Ft — e Fs '■= £f, f° r positive applied bias the tunneling from surface to tip through the 
adsorbate level is prohibited. The excitation is forbidden because all states at the tip are 
occupied up to the energy ep, thus making it impossible for an electron from the surface 
with energy e F — \eV\ — hu to tunnel to the tip. The probability of relaxing an adsorbate 
vibrations due to the inelastic scattering of tunneling electrons from the surface to the tip 
is negligibly small, because the scattered electron would need several Huu to gain enough 
energy. But this process can also be excluded, since the electron-vibration interaction on 
the adsorbate vibration is in general weak and we treat it by first-order perturbation theory. 
Thus, the transition rates can be written as, e.g. 



T^St 

l f - 



l^ff / de p° m (e - hu) fje) V \eV\ > hw 

e F -\eV\+hw 



V \eV\ < hjj 



(A2) 



,st \eV\>hw 



4 A s A* Aq 



tt A h 2 u (4A 2 + h 2 u 2 
+ tan _1 



\hu} ^tan 1 


— € m + Ep 


+ tan 1 


e m -e F + \eV\ 


A 


A 







s m — Bp — Huu + \eV\ 



+ tan 



-i 



-e m + e f — tvjj 



+ A (log [A 2 + (e m - E F f] + log [A 2 + (e m -e F + \eV\f] 
- log [A 2 + (e m -e F -hu+ \eV\) 2 } - log [A 2 + (e m - e F + huf] ) \ , (A3) 
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4 ~ Mw 



/ <fe PmO + M 



e p — \eV\—hu 



4A s A t A 2 



7T A fr 2 w (4A 2 + ft 2 u; 2 



Tjw ( tan 1 



+ tan 1 



A 



tan 



e F + \eV\ 
~A 

— e m + ep 
A 



+ tan^ 1 



-£ m + Ep + ?iUJ 

A 



+ A (log [A 2 + (e m -e F - hu) 2 ] + log [A 2 + (e m -e F + fiw + \eV\) 2 ] 

-log[A 2 + (e m -e F ) 2 ] -log[A 2 + (e m - £F + |e\/|) 2 ])| . (A4) 

The parameter Ao given in Eq. (j3J) is an important parameter in our theory, as one can 
clearly see in the Eqs. (A1-A4). In contrast to all other parameters, e.g. the broadening 
A or the energy e m , it is in general difficult to determine it from experiment or ab-initio 
calculations. However, these are only prefactors which change the absolute magnitude of 
the transition rates and thus can be easily fitted to the experiments. 
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